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We present a theory for a complete reconstruction of non-local spin correlations in ferromagnet- 
normal conductors. This quantum spin tomography is based on cross correlation measurements of 
electric currents into ferromagnetic terminals with controllable magnetization directions. For normal 
injectors, non-local spin correlations are universal and strong. The correlations are suppressed by 
spin-flip scattering and, for ferromagnetic injectors, by increasing injector polarization. 
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Spintronics utilizes the electron spin in electronics ap- 
plications and is an important subfield of condensed mat- 
ter physics. It is possible to create metallic or semicon- 
ducting hybrid ferromagnet-normal conductor systems 
smaller than the spin- flip length [l|, [2| , yet semiclassically 
large. Topics of current interest such as spin injection, 
precession, and relaxation [ll-0| , spin Hall effects [J] , cur- 
rent induced magnetization excitations |, the reciprocal 
magnetization dynamics induced spin-pumping [6|, spin 
based transistors [7j, and ferromagnet-superconductor 
heterostructures Q focus on the average non-equilibrium 
spin accumulation and dynamics. 

The correlations between injected spins in 
ferromagnet-normal conductor systems have received 
much less attention. In two-terminal junctions, current 
correlations have been investigated in few-level quantum 
dots as well as semiclassically large systems loL 11|. 
The prime targets have been noise due to spin-flip 
scattering and the super or sub poissionian nature of the 
auto correlations. 

In multiterminal junctions, current cross correlations 
allow investigations of non-local spin transport proper- 
ties. Of main interest has been the sign of the cross 
correlations, studied in quantum dots [121 ] . diffusive [13 1 
and superconducting [1J] systems and chaotic cavities 
floT ]. Moreover, in the context of entanglement of itin- 
erant spins, works on few-mode 16[ and recently also 
semiclassical 17], LL8| conductors considered non-local de- 
tection schemes with cross correlations between currents 
in non-collinear ferromagnetic terminals. 

A fundamental and important question which has not 
been addressed is if known non-local spin injection and 
detection schemes [l|, Q can be extended to identify non- 
local spin-correlations. Imagine spins injected into a nor- 
mal conductor and detected at two different spatial loca- 
tions by ferromagnetic terminals. What are the non-local 
spatial correlations between the spins? Is it possible to 
completely characterize the correlations by experimen- 
tally accessible electrical current correlations? We pro- 
vide answers to these questions for semiclassical systems: 
i) non-local spin correlations are strong, and for normal 
injectors, universal and ii) spin-correlations can be recon- 



structed by a sequence of measurements of correlations 
of currents at ferromagnetic detectors with controllable 
magnetization directions, a quantum spin tomography. 

We consider a semiclassically large, normal (metal or 
semi-) conductor connected to a normal or ferromagnetic 
injector, biased at a voltage V, and two spatially sepa- 
rated detectors, A and B, see Fig Q] Detector A (B) 
consists of a normal node coupled to grounded ferromag- 
netic terminals Al and A2 (Bl and B2) via tunnel con- 
tacts with conductances Gai and Gai {Gbi and Gb2)- 
Throughout, conductances are dimensionless and in units 
of the conductance quantum 2e 2 /h. The detectors A and 
B probe non-invasively the non-local spin correlations. 
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FIG. 1: a) A normal conductor is connected to an injector bi- 
ased at voltage V and two detector nodes A and B. The node 
A (B) is coupled to grounded ferromagnetic detector termi- 
nals Al and A2 (Bl and B2). b) Node A is connected to the 
normal conductor, as well as nodes Al and A2 via tunnel con- 
ductances Ga, Gai, and Gai, respectively. The polarizations 
Pai and Pa2 of the contacts to the ferromagnetic terminals 
are in opposite directions. 

Let us first summarize and explain our main results i) 
and ii) for the non-local correlated spin transport prop- 
erties in the device in Fig. [TJ First, combining scatter- 
ing theory and a Boltzmann-Langevin approach we de- 
rive an expression for the current correlations S A iBj — 
(2e 2 /h) Jq V dEsAiBj(E) with 

s AiB:i = 4G At G Bj ((Sf c A + V M ■ Si A )(Sf B + P Bj ■ 8f B )) f , 

(1) 

where J 3 At (P-Bj) is the polarization of the tunnel con- 
tact to terminal Ai (Bj), 5J A /b — ^Ia/b^ + ^a/b " & 
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is the fluctuating part of the 2x2 spin distribution ma- 
trix at A/B, <x = [<7 x ,<7 y ,<7 z ] is a vector of Pauli ma- 
trices, and (..)/ denotes the average over fluctuations. 
The matrix 5 Jab, with elements 8f AB — (8f A 8f B )f, 
p, q G {c, x, y, z}, is the spin correlation matrix, describ- 
ing the irreducible, or exchange, correlations between 
spins at A and B. 

We then show our result ii) : 8f AB can be reconstructed 
by a sequence of measurements of e.g. Saibi with differ- 
ent settings of P^i and Pbi- Importantly, this quantum 
spin tomography can be performed for arbitrary (finite) 
magnitudes of the polarizations |Pai| and |Psi| and 
spin-flip scattering in the conductor. Moreover, global 
spin symmetries limit the number of finite elements of 
SfAB-, allowing for a simplified quantum spin tomogra- 
phy with fewer cross correlation measurements. 

For a normal injector we derive a generic expression 
for 5Jab, with nonzero elements 

5 Sab = Sf AB /2, Sf% B = 8ff B = 8ff B = jSf AB /2 7 

(2) 

where 8f AB is the equal-spin correlator and 7 quantifies 
the spin coherence in the conductor. 7 = 1 for a coherent 
system, i.e. no spin-flip scattering, and 7 = for a sys- 
tem with strong spin-flip relaxation. For a ferromagnetic 
injector, the correlations depend on the properties of the 
conductor, as shown below. 

Inserting Eq. ([2]) into (JT]) gives a cross correlator 



({k} = {c,x,y,z}) 



s AiBj = 2G Ai G B jSf AB [l + 7 P A i ■ P B j], 



(3) 



depending on the relative orientation of the polarizations 
P A i and Psj- This together with Eq. ([2]) demonstrate 
our counter- intuitive result i): any conductor with a nor- 
mal injector displays strong and universal non-local spin- 
correlations. We note that for the current cross correla- 
tor, similar results have been obtained in particular ge- 
ometries 16- IH with no spin- flip scattering, 7=1. 

We now describe the quantum spin tomography, start- 
ing for clarity with the known properties [IpJ of the 
average spin distribution matrix in node A, f A = 
f A l + f A ■ <r, where the real polarization vector = 
[f A ,f A ,f A ] with |f A | < 1. The average current is I A i = 



(e//i) Jq V dEi A1 (E) with [l9| 
i A \ = 2Gai [/a 



ai 



(4) 



For the quantum spin tomography, we transform the or 



bital scheme developed in Ref. [21( to the spin degree 
of freedom and extend it to account for arbitrary de- 
tector polarization. Formally, to determine f A ,f A four 
independent measurements of the current are needed. 
The theoretically most convenient set {I AX }, k = 1 — 4 



has the polarizations P A l/P A i = [0, 0, 1], P A l/P A i — 
[0,0,-l],P^/^i = [1,0,0] and Pf x /P A i = [0,1,0], 
where Pai = |Pai|, but other settings are also fea- 
sible. The expression in Eq. Q then allows writing 



(2) 
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(5) 



Knowing the polarization P A i and the conductance Gai 
from independent measurements, the spin-distribution 
matrix j a is fully reconstructed by current measure- 
ments. Importantly, for a normal injector, only f A is 
non-zero. For a ferromagnetic injector, when the spin 
quantization axis along the direction of polarization, only 
f A and f A are non-zero. 

We then turn to the spin correlation matrix <5/ab, with 
the 16 real elements §f v AB - This implies that we need 
16 independent cross correlator measurements to deter- 
mine all elements $f AB and reconstruct 6 Jab- From Eq. 
^ we obtain the formal relation between the coefficients 
Sf AB and the cross correlators 



Sf pq 



AB 



8G A iG B iVe 3 /h ^ 



E n pk n ql q(k,l) 
{ ^A1^B1°A1B1 



(6) 



where S A { B1 is the cross correlator with the detector ter- 
minal setting k at Al and I at Bl. Here Qbi is obtained 
from Qai by changing P AX to Pbi- 



For a normal injector, the requirement 22j of invari- 
ance of 5 Jab under any global spin rotation means that 
there is only four non-zero elements 8f AB and Sf AB = 
^ fjU3 ~ ^f A B' F° r a ferromagnetic injector (defining 
the spin quantization axis) invariance of 5f AB under 
the global rotation {|t),|4-}} -> {e l0 |t), e"* ll)} yields 
[23I ] six non-zero elements S f c AB , 8 f AB , S f AB , S f AB and 

SfTB = Sf y A y B- 

From Eqs. §5§ and © the detector polarization set- 
tings necessary to determine the non-zero components 
of /a and 8 Jab are found: For a normal injector, only 
collinear polarizations at A and B are needed for both j a 
and 8 fAB- For a ferromagnetic injector, for Ja the de- 
tector polarizations in addition have to be collinear with 
the injector one. However, for 8 Jab non-collinear polar- 
izations at A and B are necessary, e.g. both along the x 
and z axis, since 8f AB ^ Sf AB = Sf A v B . Importantly, for 
an unknown direction of the injector polarization or two 
(or more) non-collinear ferromagnetic injectors, the full 
tomographic scheme with detector polarizations along all 
three axes x,y and z are required. 

We will now detail our calculations, assumptions, and 
approximations. In addition to the information given 
above, the normal conductor in Fig. [T] is connected to 
detector nodes A and B via tunnel barriers with conduc- 
tances Ga and Gb- The two ferromagnetic terminals Al 
and A2 (Bl and B2) have opposite directions of polariza- 
tion. We assume the limit of low temperature kT <C eV. 
All conductances are much larger than unity. 
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It is assumed that the normal conductor consists of 
diffusive and/or chaotic parts, allowing a semiclassical 
treatment of the orbital properties. In contrast, spin is 
treated fully quantum mechanically. Furthermore, scat- 
tering is clastic. Following the magnetoelectronic circuit 
theory of Ref. [19|, we discretize the system into nodes 
connected via tunnel barriers, see Fig. [TJ Each node v, 
spatially much smaller than the spin-flip length, is char- 
acterized by a 2 x 2 distribution matrix with an average, 
f v , and a fluctuating, 8f v , part. To ensure that the detec- 
tors do not influence the spin-properties of the system, we 
require i) Ga < Gai+Ga2 and Gb < Gbi + Gb2 so that 
an electron entering e.g. node A from the conductor is 
emitted into Al or A2 and do not return to the conductor 
and ii) G^iPai = ~G A 2Pa2 and GbiPbi = -G B 2Pb2, 
which ensures that no spin polarization is induced into 
the conductor from the ferromagnetic terminals, i.e. the 
measured spin signal arises from the conductor exclu- 
sively and not from the detector circuits. 

Deriving Eqs. (|4|) and ([1]), we first review [3] the 
spin information present in the av erag e spectral current 
%ai (E) . In the scattering approach [201 ] , with no particles 
incident from terminal Al in the bias window (0 < E < 
eV), the spectral current is 



lAl 



L Al, 



°Al.n°Al, 



(7) 



where b A \ creates an electron on the ferromagnetic side 
in the contact between Al and A, in conduction mode n 
propagating into Al and the energy-dependence is sup- 
pressed. The spin quantization axis a = j" , 4- is along the 
direction of Pax- The creation operators b A \ n are re- 
lated to the operators b Am for electrons on the normal 
conductor side, emitted from node A towards Al, via 
the spin-dependent transmission matrix of the normal- 
ferromagnetic interface tAi with elements t a A \ nm . Fol- 
lowing Ref. [jjj], we make the semiclassical approxima- 
tion that the spin distribution matrix in node A is in- 
dependent on mode index, i.e. (6°^ 
giving 



J An^Am) 



J A u nm, 



iAi =J2 f A?fl T = GUitr {(1 + P A i ■ ar)h} ■ (8) 



Here [Tj Tax = E„,„(^i,„m) T (tAi,r im )* = G A1 {1 



where n T Blm = b B \ m b T Bl m and b T B \ m creates an out- 
going electron on the ferromagnetic side, in conduction 
mode m in Bl with spin quantization axis along the di- 
rection of the magnetization ng. Disregarding terms of 
second order in Ga/{Gai + Ga2) or Gb/{Gbi + Gb2), the 
operators b A \ n and b T B \ m are expressed in terms of the 

operators b A l and b B J and the scattering amplitudes of 
the respective normal-ferromagnetic interfaces. Making 
the semiclassical approximation that the non-local irre- 
ducible correlator (b A ib^b B \b B \}-(b a A lb^ Am )(b B \b B \) = 

SfAB ,TT Snmfah we arrive at 



SA1B1 — 'Al 'Bl JAB —^A\^B\ 

era' tt' 

x tr { [(i + P A i ■ &) ® (1 + Pbi ■ a-)] 5f AB ] (10) 

where Tbi is obtained by changing all indices A to B in 
Tax and ® is the tensor product. Here we work in the 
basis {|ot)} = [| ft), | tl), I It), I II)], i-e- the matrix el- 
ements (5fAB)\<jT),\<j'T') — Sf AB ' TT ■ As is shown below, 
the 4x4 spin correlation matrix S/ab — {SJa <& Sfs) /, 
provides a semiclassical interpretation of SfAB- This 
means that Eq. (jTUJ) directly gives Eq. Jl]). Moreover, 
the expressions for the other correlators SAiBj can sim- 
ilarly be given in terms of SfAB- This shows that 8 Jab 
contains all information about non-local spin- correlations 
that can be obtained from cross correlations. 

To further investigate the properties of f a , fB and 
SfAB we now turn to the spin-dependent Boltzmann- 
Langevin approach of Ref. llj. The average part of 



the distribution matrix /„ at node v is determined from 
the condition of conservation of matrix currents into the 



node , 



0, iyn 



The 2x2 matrix current 



between a normal node v and a ferromagnetic or normal 
node [x is i vyL = (G„ M /2) |(i + P M • &), (/„ - /„) J with 

{.,.} the anti-commutator, G v>Jl the tunnel conductance 
between the nodes and P n the polarization vector of node 
(i (P n = for a normal node). The distribution matrices 
for normal and ferromagnetic terminal nodes are 1 for 
biased terminals and for grounded. This allow us to 
calculate the distribution matrices of all nodes. 

For the fluctuating part of the distribution matrix, we 
first note that the total fluctuations of the matrix current 



(t 



Al 
Al 



&) where the elements of the' 2 x 2 matrix are fl owing between two nodes v and /i is a sum of the 



nm jar 



4-CTT 

b Al,nm- 



Eq. 



directly gives Eq. ((4]). 



Similar relations hold for the average currents into A2, 
Bl and B2. 

We then turn to the low frequency correlations between 
electrical currents in e.g. terminals Al and Bl, Saibi — 
J dt{AI A i{0)AI B1 (t)). Scattering theory [2(| gives 



SA1B1 — 



[( n Al.n n Bl,m) ~ ( n Al,n)( n Bl,m)] ( 9 ) 



bare fluctuations 5i vli and 5T VU due to the fluctuating 
distribution matrices. For v normal and fi normal or fer- 
romagnetic 5i ytl = (G„ M /2) j(i + P n ■ &), (<5/„ - j. 
The requirement of matrix current fluctuation conser- 
vation Ai va = then gives 5f v in terms of 8i ua . 

The bare fluctuations di ufl at different contacts are 
uncorrelated while for v, /i normal {8i v ^ ® Si^ u ) / = 

(G^/2) [/„ ® (i - /„) + U ® (i - M] w + h - c - where 



4 



h.c. denote hermitian conjugate and the permutation ma- 
trix W has nonzero elements Wxi = W23 = W32 = 
W44 = 1. For v normal and fj, ferromagnetic we have 
(fiu„ ® f = (G^/2){([1 + P„ ■ ® (1 - f v ) + 

([i +P M • <r][l - /„]) (g) /„}VF + /i.e.. Here we used that 
ferromagnetic (i.e. terminal) distributions do not fluctu- 
ate. From these relations any electrical current correlator 
(Ai^Ai^) f, with Ai v — tr[Az„], can be obtained. 

Spin flip scattering is taken into account on the level 
of the relaxation time approximation. This amounts to 
coupling each node n to a spin-flip node <pu with a tun- 
nel contact with conductance G vv oc 1/ T vv , with r v „ the 
spin-flip time of the node, and requiring conservation of 
electrical current and current fluctuations into the spin- 
flip node. Here we give the universal results of the cal- 
culation, i.e. we consider an arbitrary normal conductor 
with any amount of (spatially dependent) spin-flip scat- 
tering, the details of the calculations are given elsewhere. 
First, by comparing the obtained expression for the spec- 
tral cross correlators SAiBj = {AiAiAiBi) / with Eq. (fT0|) 
we conclude that 5 Jab — (SJa ® Sfs)f, discussed above. 
Second, for normal injectors, we find the generic form 



Sf AB = (Sf%/2) (l- 7 )i®i + 27W 



(11) 



This is just the result in Eq. ([2]). 

Further insight is obtained by calculating the prop- 
erties of the simplest possible conductor, a single node 
0. For a normal injector we find the distribution func- 
tion at e.g. A as Ja = Ga/(Gai + G^/l, with 
/ = G/(G + Ga + Gb) the distribution function of the 
conductor node and G the injector-conductor node con- 
ductance. This is independent on spin-flip scattering. 
For the spin-correlation matrix we get the result in Eq. 
m with Sf AB = f 3 G A GB/[G(GAi + GA2)(GBi + G B 2)} 
and 7 = [1 + t/t^- 1 with t/t v = G V I(G A + G B +G) 
the ratio of spin-flip and dwell times in the central node. 

For a ferromagnetic injector with polarization P/ the 
spin distribution matrix at e.g. A has two non-zero com- 
ponents (f A ) n = f A and (f A ) U = f A with = 
f[(l ± Pr)(l =F P t f) + t v /t]/[(1 - Pff) + t v /t] with 
Pi = |P/|. For the spin-correlation matrix, the full 
expression, including spin-flip scattering, becomes very 
lengthy and we only present the result for 7=1. This 
is 6/%% = SfX% = Sf AB (c + + c_)/2, Sff B = 5f%% = 
SSf AB (c + - c_)/2 and Sf%% = 
c± = (l±Pr) 2 /(l±-Pr/) 3 and c 

Inserting this into Eq. (jlOl) we get the cross correlator 



AB = tf°ABC0/2 With 
(1-P|)/[I-(P / /) 2 ]. 



Saibi = 2GAiG A iSf A B(ci + co[l + P 



Al 



(12) 



with CI = (1 + P z Al P z Bl )(c + + c_ - 2c )/2 + (P z Al + 
P B1 )(c + — c_)/2. This clearly demonstrates that while 
a ferromagnetic injector leads to a polarization of the 
conductor, it suppresses the spin-correlations. 

In conclusion we have presented a scheme for quan- 
tum state tomography of non-local spin correlations in 



normal-ferromagnetic conductors. Non-local correlations 
are generically strong but suppressed by spin-flip scatter- 
ing and ferromagnetic injectors. 
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